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We report the Standard Model evaluation of the indirect CP violation parameter ek using inputs 
determined from lattice QCD: the kaon bag parameter Bk, ^o, |14is| from the K 13 and decays, 
and |14f,| from the axial current form factor for the exclusive decay B —>■ D*£i' at zero-recoil. 

The theoretical expression for ek is thoroughly reviewed to give an estimate of the size of the 
neglected corrections, including long distance effects. The Wolfenstein parametrization (IJkbl, A, p, p) 
is adopted for CKM matrix elements which enter through the short distance contribution of the box 
diagrams. For the central value, we take the Unitarity Triangle apex (p, p) from the angle-only ht 
of the UTfit collaboration and use Vus as an independent input to fix A. We find that the Standard 
Model prediction of ek with exclusive Vcb (lattice QCD results) is lower than the experimental 
value by 3.4a. However, with inclusive Vcb (results of the heavy quark expansion), there is no gap 
between the Standard Model prediction of ek and its experimental value. For the calculation of 
Ek, we perform the renormalization group running to obtain rjee at next-to-next-to-leading-order; 
we find 7 ?™^° = 1.72(27). 
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I. INTRODUCTION 

CP violation in nature was first discovered in an ex¬ 
periment with the neutral kaon system in 1964 [1] . There 
are two kinds of CP violation in the neutral kaon system: 
one is the indirect CP violation due to CP-asymmetric 
impurity in the kaon eigenstates in nature, and the other 
is the direct CP violation due to the CP violating nature 
of the weak interaction [HE]. CP violating observables 
are prime candidates in searches for physics beyond the 
Standard Model. Experimentally, CP violation in the 
neutral kaon system is known more precisely than in any 
other physical system. Here, we focus on the indirect CP 
violation in neutral kaons. 

Indirect CP violation in the neutral kaon system is 
parametrized by Ek 

^ -> 7r7r(/ = 0)) 

^ A{Ks -T 7r7r(/ = 0)) ’ 

where and Ks are the neutral kaon states in nature, 
and 1 = 0 represents the isospin of the final two-pion 
state. In experiment |4], 

£k = (2.228 ± 0.011) X 10"^ X , 

4 = 43.52 ± 0.05° . (2) 

Here, the Ek value represents an « 0.2% impurity of the 
CP even eigenstate in the Kl state, which contains 99.8% 
of the CP odd eigenstate. 

We can also calculate Ek directly from the Standard 
Model (SM). In the SM, CP violation comes solely from a 
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single phase in the CKM matrix elements EUSj. TheSM 
allows the mixing of neutral kaons K'^{ds) and K {sd) 
through loop processes, and describes contributions to 
the mass splitting AMk and Ek- Hence, we can test the 
SM through CP violation by comparing the experimental 
and theoretical values of Ek- 

In the SM, the master formula for ek is 


Ek =e^^V2sm9^CeBKXsD + ^ + Cld) 

+ 0{u;E') + 0{^or2/ri), (3) 


where Cg is a well-known coupling given in Eq. (68), and 


XsD is the short distance contribution from the box dia¬ 


grams given in Eq. (61). Here, the major contribution to 


Ek comes from the Bk term, and the minor contribution 
of about 5% comes from the ^0 term. The remaining 
contribution of ^ld is about 2% coming from the long 
distance effect on ek [71 [Hj. In Section [h{ we re-derive 
the leading contribution given in Refs. |9tillj. We also 
explicitly derive higher order corrections, including the 
long distance contribution, in this paper. A similar for¬ 
mula without the long distance correction ^ld and higher 
order terms appears in Refs. pnuj. 

In order to calculate Ek directly from the SM, we use 
input parameters obtained from lattice QCD and exper¬ 
iments. In particular, Bk and Vet have dominated the 
statistical and systematic uncertainty in the SM evalua¬ 
tion of Ek for a long time. 

During the past decade, lattice QCD has made signif¬ 
icant progress in calculating Bk so that its error is re¬ 
duced dramatically, down to the « 1.3% level at present. 
This result is available from the Flavour Lattice Aver¬ 
aging Group (FLAG) [Hj. It is obtained by taking an 
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average of the Bk results from a number of lattice QCD 
groups [EHn]. We calculate Ek using two different input 
values of 13k- one is the FLAG result [12], and the other 
is the most updated result from the SWME collaboration 

m- 

It is also noteworthy that the lattice calculation of the 
amplitude ImA 2 related to the decay K —>■ = 2) 

m makes it possible to determine more precisely. 

Another important input parameter to Ek is Vcb- 
There are two independent methods to determine Vet- 
one is the exclusive method j^D], and the other is the 
inclusive method |2T1122| . In the exclusive method |20| , 
one uses lattice QCD to calculate semileptonic form fac¬ 
tors for the decays B —>■ In the inclusive method 

ED [52], one performs analysis on i? —> decay pro¬ 

cesses using the heavy quark expansion [23] ■ Here, we use 
both the exclusive and inclusive V^b to determine Ek^ and 
we compare the results with each other and experiment. 

We use the Wolfenstein parametrization for the CKM 
matrix E3], truncating the series at O(A^) « 10“^. Here, 
we use three different choices of Wolfenstein parameters: 
(1) A, p, and fj from the global fit of the CKMfitter col¬ 
laboration [25l|26], (2) A, p, and fj from the global fit of 
the UTfit collaboration ED EH], and (3) p and rj from an 
angle-only fit (AOF) from the UTfit collaboration [55] . 
with an independent input for A directly from Vus [I]- 
In all the cases, we take Vcb instead of the Wolfenstein 
parameter A from the unitarity triangle (UT) analysis. 
We emphasize that the AOF does not use Ek, Bk, and 
Vcb to determine the UT apex p and fj. Hence, it pro¬ 
vides a self-consistent way to test the validity of the SM 
with Ek, using the lattice results for Bk and Vcb with no 
correlation between {I3K,Vcb) and {p,fj). 

To estimate the effect of correlations in lattice input 
parameters, we note that Vcb dominates the error in ek, 
and the FLAG 13k [H] is dominated by the BMW collab¬ 
oration result m- We assume that there is no correlation 
between the BMW Bk and the exclusive Vcb from the 
FNAL/MILG form factor [5D], because their gauge en¬ 
sembles are independent. Hence, we assume that the cor¬ 
relation between the FLAG 13k and the FNAL/MILG Vcb 
are negligibly small. However, when we use the SWME 
33K |18j . there must be an inevitable correlation with 
the ENAL/MILC result for exclusive Vcb- In this case, 
we consider -1-50% correlation and —50% anti-correlation 
between the SWME Bk and the exclusive V/b to esti¬ 
mate the systematic error due to the correlation between 
them. The RBG/UKQCD collaboration calculated us¬ 
ing domain-wall fermions, which is also completely inde¬ 
pendent. Hence, we assume that is uncorrelated with 
the other lattice inputs 13k and Vcb- 

When we determine the value of Ek, we take into ac¬ 
count the correlation between the SWME I3k and the 
FNAL/MILG Vcb and assume that the other input pa¬ 
rameters are uncorrelated. We use the Monte Carlo 
method to calculate the Ek distribution from the SM. 
The results are cross-checked using the standard error 
propagation method. 


In Section [TTJ we review neutral kaon mixing and de¬ 
rive the master formula for ek from the SM. Here, we 
give an estimate for the size of truncated small correc¬ 
tions. In Section we explain each input parameter 
in detail. Here, we also explain details on how we pop¬ 
ulate input distributions using the Monte Carlo method 
and how we determine errors on Ek considering differ¬ 
ent input combinations and correlations among them. In 
Section |TV[ we present the results for ek obtained using 
various combinations of input parameters. In Section jV] 
we conclude. 


II. REVIEW OF EK 

A. Effective Hamiltonian 

Let us first review the theoretical formalism of neutral 
kaon mixing in the SM [5D] . Let us consider a state that is 

initially (at t = 0) a superposition of K°{ds) and K^{sd): 

1^(0)) =a(0)|A°) +6(0)1 a"). (4) 

This state will evolve in time, and part of it will decay 
into final states {/i, / 2 , ...} as follows. 

m)) = a{t)\K^)+b{t)\K°) 

+ Cl(Q|/l) +C2(t)|/2) H-• (5) 

If we are interested in calculating only the values of a{t) 
and b{t), but not the values of Ci(t), and if the time t 
is much larger than the typical strong interaction scale, 
then we can use the simplified formalism in Ref. mm- 
In this formalism, the time evolution is described by a 
2x2 effective Hamiltonian HeS that is not Hermitian, 
which allows the neutral kaons to oscillate and to decay. 

The neutral kaon system forms a two dimensional 
subspace of the Hilbert space of the total Hamiltonian 
H = Hq -f- iLw Hq is the strong interaction Hamiltonian 
which defines the full Hilbert space. Decays into different 
strong eigenstates are mediated by the weak interaction 
Hamiltonian which is treated as a perturbation. 

In the 2-dimensional subspace, the time evolution of 
the neutral kaon state vector can be described by the 
effective Hamiltonian iLeff, 

tj^\K{t)) = Hcs\Kit)) - (6) 

The effective Hamiltonian consists of two Hermitian op¬ 
erators M and F, 

Hc« = M-i^. (7) 

The dispersive part M defines masses of the neutral kaon 
states, which correspond to the kaon eigenstates in na¬ 
ture, and the absorptive part F defines decay widths of 
the mass eigenstates in the presence of the weak interac¬ 
tion iLw The effective Hamiltonian itself, however, is not 
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Hermitian. It is a necessary consequence to take into ac¬ 
count kaon decay amplitudes that have final strong eigen¬ 
states which do not belong to the neutral kaon subspace, 
as one can see in Eq. (§. 

The decay processes can be systematically described 
by the perturbative corrections to the effective Hamilto¬ 
nian of the subspace [33]. In the second order in or 
equivalently second order in the Fermi coupling constant 
Gf, the results are, as shown by the famous Wigner- 
Weisskopf formula [^ |32| , 


{a\H^\C){C\H^\(i) ^ 

^ Ec- rriKO 

=2fY,{c^\H^\C){C\H^W{Ec - mo), (9) 

c 

where mxo is the mass of the neutral kaons if® and if°, 
V denotes the principal value, \C) is an intermediate 
state with energy Ec which belongs to the full Hilbert 
space, and the summation over C includes integration 
over the continuous quantum numbers. Here, we ignore 
a tiny experimental mass difference between if® and if , 
since we assume CPT invariance throughout this paper. 
Hence, the masses of a particle and its anti-particle are 
the same. 

The leading correction to the off-diagonal components 
of Map comes from the four-quark AS = 2 operator of 
dimension 6. It is built from a product of two weak 
current-current interactions by integrating out VT-bosons 
and heavy quarks in the box loop diagrams. This is a 
short distance contribution, and it is the leading effect 
which is responsible for the indirect CP violation in neu¬ 
tral kaon mixing. This short distance effect is explained 
in Section III Cl in detail. If there exists a fundamental 
AS” = 2 interaction, the so-called superweak interaction 
ifsw, which is absent in the SM, it also contributes to 
the off-diagonal components M^p [34]. Neutral kaons 
could decay into an intermediate state \C) as a result of 
AS” = I transitions. The parts which involve these inter¬ 
mediate states \C) in Eq. ([^ and Eq. ([^ constitute the 
long distance contributions. 

The time independence of the effective Hamiltonian is 
a consequence of the Wigner-Weisskopf approximation, 
which takes the interaction time to infinity and turns the 
interaction adiabatically off [33l [35] . The well-known ex¬ 
ponential decay law follows from this approximation. So 
a deviation from the conventional exponential decay gives 
an estimate of the accuracy of the Wigner-Weisskopf ap¬ 
proximation. These corrections to the exponential de¬ 
cay, with present and foreseeable experimental precision 
[551155] , are far beyond the precision that we pursue here 
for the value of Ek in the SM. Hence, we neglect these 
corrections in this paper. 

Before considering explicit calculation of the matrix 
elements on the right hand side of Eq. ([^ and Eq. 
we focus on their parametrization. From the Hermiticity 


of M and P, each of them is parametrized with 4 real 
variables 


f Ml im' + 6 m'\ 

\-im' + 5m' M 2 ) ’ 
f Pi *7' -I- 5^' 

\^—iy'-I-r2 



Further simplification 


^ m ' — 0 7 ^7' — 0 


( 10 ) 

( 11 ) 


( 12 ) 


follows from CPT invariance, {CPT)0{CPT)~^ = O, 
where O = M, P, and with a specific basis made of the 
CP eigenstates: {\Ki), \K 2 )} [55] . 

Assuming the strong interaction has CP symmetry, the 
neutral kaon subspace can be spanned by the CP even 
\Ki) and odd \K 2 ) eigenstates 


\Ki)=y{\K^)-\K^)), 

\K2)=y{\K^) + \K^)). (13) 

We adopt a phase convention [50] 

CP\K°) = -|Z°) , (14) 


and for time reversal T 

T\K°) = -|iP®), r|]p°) = . (15) 


Here, note that the incoming state becomes an outgoing 
state under time reversal and vice versa. Then 


\Ki) =CPT\Ki) = -\Ki), 

IZ 2 ) =CPT\K2) = \K2) . (16) 

Then we can verify the constraints in Eq. (HI), 

{Ki\M\K2) ={K2\{CPT)M\CPT)-^\Ki) 

= -{K2\M\Ki). (17) 

The same relation also holds for P. Here, the Hermi¬ 
tian conjugate arises from the anti-unitarity of the time 
reversal symmetry. 


B. Ek and i 


The physical states Ks and Kl are approximately CP 
even and odd, respectively. In other words, the physical 
eigenstates of the effective Hamiltonian TPeff in Eq- 0 in- 
elude a tiny impurity (« 10“®) of the opposite CP eigen¬ 
state defined in Eq. (13). The physical eigenstates can 


be written with small mixing parameters is and el, 


\Ks) = 


1 


\/l + l^sl' 


,{\Ki)+es\K2)). 
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\Kl) = , i\K^) + e^m ). ( 18 ) 

Their eigenvalues are 

\s = A — AA, Al = a + AA, (19) 

where 

A = ^|(Mi + M2)-^(ri + r2)| , (20) 

AA = ^y (am + 1 at)' + 4(m' - ^ 7 ')' , (21) 

and 


To complete the connection between e and ek, we 
need to consider kaon decay amplitudes m- Define the 
isospin amplitude Aj and phases and Sj by 

AiK° 7 r 7 r(/)) = . (29) 

Then, in our phase convention, which is one of the most 
popular conventions |55] . 

A(K° 7 r 7 r(/)) = - A*jA^‘ = -\Ai\e-^^‘A^^ , (30) 

and 

A{Ki 7r7r(/)) = V2ReAiA^‘ , 

A{K 2 7 r 7 r(/)) = iV2\mAiA^‘ , (31) 


AM = M 2 - Ml, AT = Ti - r2 . (22) 

Eliminating the eigenvalues Xs,l from the system of 
eigenvalue equations 

(^Mi - ^Ei - As^ + es (im' + ^ 7 ') = 0 , 
es(^M 2 - ^r 2 - As) - (im' + ^ 7 ') = 0 , 

(m 2 - ^r 2 - Al) - El (im' + ^ 7 ') = 0 , 

SL (Mi - "-Ti - Xl) + (im' + ^ 7 ') = 0 (23) 

leads to the condition 

(4,l + 1) (im' + ^ 7 ) - £s,L (am + ^ Ar) = 0 (24) 

that the mixing parameters have to satisfy. The 

quadratic equation in Eq. ( p4| has two solutions. One 
of them is very small (« 10“^) and the other is very 
large (« lO'*'^). Hence, it is obvious that the two mixing 
parameters are equal, since we assume that the mixing 
impurity is in the level of « 10“^. Hence, we will use the 
mixing parameter i 


£ = £S = £l- 


(25) 


Since we know that |e| 
as follows. 


10 we can rewrite Eq. (24) 


£ — ^(o)(l + £^) I (26) 


where 


where the phase Sj is equal to the S'-wave scattering 
phase shift of the final two-pion state by the strong inter¬ 
action, and the subscript / represents the isospin of the 
final state. Assuming isospin symmetry, this comes from 
Watson’s theorem [inET]. Watson’s theorem is based on 
time reversal symmetry implicitly. Because the final state 
scattering only involves Hq, application of Watson’s the¬ 
orem concerns the time reversal symmetry of the strong 
interaction. It is equivalent to the CP symmetry, if we 
assume CPT invariance. Here, note that represents 
the effect of the violation of Watson’s theorem. 

If the weak Hamiltonian 77^ respected CP symmetry, 
which is equivalent to time reversal symmetry under CPT 
invariance, then Watson’s theorem must hold to guaran¬ 
tee that A] must be real in this case m- However, we 
know that i7w breaks CP symmetry through the exis¬ 
tence of a single phase in the CKM matrix, and so it also 
violates time reversal symmetry. As a consequence, Wat¬ 
son’s theorem is violated, and so Aj becomes complex, 
which generates the phase ^7 ^ 0, in general. Hence, the 
weak phases parametrize the direct CP violation in 
the weak interaction with a non-zero phase difference, 
Im(A 2 /Ao), which is independent of phase convention 

m- ^ 

Now, let us focus on 7 ' and AP in Eq. ( 27). We 
will address m' and AM later in Section III Cl and Sec¬ 


tion [H^ Let us divide both numerator and denominator 


of Eq. (271 by AM; we obtain 


£(0) = e*®sin 


>( 


m . „ 7 \ 

- 7 -— — i cot 0-— ) 

AM AT/ 


sin 6 » ( - i^o cot d) - 1 - O (we') 

+ O(eor 2 /ri), (32) 


£( 0 ) = 


i(m' 


-V 
2 ' 


AM- 


lAP 


(27) 


Then, Eq. (|26|) can be solved iteratively near the leading 

(28) 


order solution e(o), 

£ = £(o) + £(o) + 


3 , 9~5 

^£(0) 


5£[o) 


where 


tand = 


e = e 


2AM 
AT ’ 

AS2-So)i^ _ ImAc 


(33) 


1/2 \ HeA 2 HeAo 


(34) 
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and 


w = 


ReA2 

Rg-Aq 


ImAo 

RgA.q 

ImA.2 

ReR? 


= tan(Co) = Co + O(Co'), 
= tan(6) = 6 + 0(C|) • 


(35) 

(36) 

(37) 


Here, we use the small angle approximation for the weak 
phases and ^ 2 - 


When we derive Eq. (32), we apply the following ap¬ 
proximation: 

^ = *^0 + O{ios') + o(6r2/ri). (38) 

It is obtained from the fact that the neutral kaon decay 
amplitudes are dominated by the 1 = 0 two-pion final 
state. First, we can express it as follows. 


ij 

AT 


= ^(i + (r 2 /ri) + (r 2 /ri )2 + ...) . ( 39 ) 

i 1 


Since we know that r 2 /ri 10 we can introduce the 
first approximation as follows, 


Z7 ij 

AT ^ ' 


(40) 


Using the Wigner-Weisskopf formula in Eq. ([^, we can 
re-express the right-hand side as follows. 


n' EciKi\H^\C){C\H^\K2)S{Ec - m^o) 

El J2c{Ki\H^\C){C\H^\Ki)SiEc - tuko) ' 


(41) 


Here, it is obvious that the denominator is completely 
dominated by the two-pion states. In the case of the 
numerator, there are contributions from two-pion states, 
three-pion states, and so on. Here, we assume that the 
two-pion contribution is dominant and we may neglect 
the rest, which includes the (semi-)leptonic decay modes. 
For example, in the case of the three-pion state, the 
branching ratio between the two-pion decay and three- 
pion decay of is about 3.5 x 10“^, and that for the 
Kl is about 113 |1]. Therefore, the three-pion decay 
mode is suppressed by a factor of about 6.3 x 10“^ com¬ 
pared to the two-pion mode. Similarly, we also assume 
that the semi-leptonic and leptonic decay modes are so 
suppressed that we may neglect them in the numerator, 
as in Refs. [lOlIM] . 

Therefore, as a very good approximation, we assume 


that the summation in C in Eq. (41) is completely dom¬ 
inated by the two-pion states in both the numerator and 
the denominator as follows. 

A{Ki ^C) = (C|Hw|iEi) 

= i5c,777r(o)(7r7r(0)|i7w|7Ei) + 


— i^c,7r7r(o) V^(ReAo)e*^“ -I-dc^7r7r(2)'\/2(ReA.2)e 


Similarly, 


182 


(42) 


A{K2 ^C) = {C\H^\K2) 

= l5c,777r(0)(7r7r(0)|i7w|7E2) -f (5c'.7r77(2) (7’‘7r(2)|iJw|7^2) 


= i5c.7i-7r(o)*v^(ImAo)e*^“ -I-dc'.7r7r(2)*v^(ImH2)e 


iSn 


(43) 


Using Eqs. (42), (43), and (41), we can obtain the fol¬ 
lowing result: 

i'y' 2z(ReHo)(IniAo) + 2i(ReA2)(ImH2) -)-••• 

id ^ “ 


2(ReHo)2 -t 2(ReH2)2 -t • • • 

ImAg f ReA2 7 f ImA.2 ImHg 

ReHn 


1 ReH2 


-f ■ 


Re Aq V Re Aq ^ 

= + ■■■ . (44) 

Here, we know that ^0 ~ 10“"^ and y/2uje' ~ 10“"^. Hence 


we may safely neglect the ioe' term in Eq. (44) within the 
precision that we pursue in this paper. This leads to the 
approximation in Eq. (38). 

In terms of isospin amplitudes, £k in Eq. ([^ can be 
written 


£k = 


HmHo -I- eReHg e -f 


ReHg -I- iHmHo 1 + 

= (e -I- *^o)(l ~ *£^0 + • • ■) 

= e -f ~ + ■ ■ ■ 

= £(0) + + £(0) ~ *£( 0)^0 + • • • 

= £(o) -f + C>(^(0)) • 


(45) 


Finally, using Eq. (32) and Eq. (K5|, we obtain 


£k = e*®sin6»( 


AM 
o(6r2/ri). 


Co) +0(c^£') 


(46) 


Here, we keep only the first two terms from Eq. (45). The 


size of those corrections that we neglect in this paper is 
much smaller than the experimental precision of Ski as 
one can see in Eq. (§. 

C. Short Distance Contribution 

The matrix element m' can be calculated from the 
Wigner-Weisskopf formula given in Eq. (|^. A short dis¬ 
tance contribution rngQ to m! is 

(6)1 


2mKO ■ im'sD = (7511177^1^1^:2) 
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s W d 

■< fWWW\A/f ^ 


U,C,t ' ' ' ' u,c,t 


> •WWWV\/> > 

d W s 


s u,c,t d 



(a) 

FIG. 1. Box diagrams for the mixing. 


=-zlm(K°\n^i^\K^) . (47) 

The factor 2m^o comes from the normalization condition 
for the external kaon states. 

In the SM, the Hamiltonian density represents 

the leading short distance term in the AS = 2 effective 
weak Hamiltonian, which is constructed from the box 
diagrams. For a scale below charm quark threshold fj, < 
He « 0{mc), 

■Hsd = ^^M^[>^lr]ecSo{xc) + XtmSoixt) 

+ 2XcXtVctSo{xc,Xt)]b{H)OLL^'^{H) + h.c. . 

(48) 


Here, the dimension-6 local four fermion operator which 
comes from the well-known box diagrams in Fig. is 

OlI^'^{h) = S7^(1- 75 )^ 57 '"(l-75)d- (49) 


By integrating out the heavy degrees of freedom in the 
loops of the box diagrams, we obtain the Inami-Lim func¬ 
tions [39] as follows, 


So{xi) 

SQ{xi,Xj) 


= Xi 


1 9 3 3a;f In Xi 

4 4(1 - Xi) 2(1 - Xi)2 (1-Xi)3_ 


XiX,- 


Vi tVj 


_4 2(1-Xi) 4(l-x*)2 

- (* ^ j) - 


iXiXj 


4(1-Xi)(l-Xj) 


In Xi 


(50) 


where, i = c,t, Xi = rrii = mi{mi) is the scale 

invariant MS quark mass m, and Mw is the IH-boson 
pole mass. The u-quark contribution is rearranged into 
c and t terms by imposing a unitarity condition, 

-l- Ac -f At = 0, (51) 

X, = V*V,d, 


and then the effective Hamiltonian Hgp is re-expressed 
with c and t terms. In Eq. (50), an approximation 
ml/M^ = 0 is used. Each pair of vertices for FF-boson 
interchange gives the products of the CKM matrix ele¬ 
ments Ai = V*.Vid. 


Besides a zeroth order QCD effect dealt with by 
the Inami-Lim functions Sq, r]ij with i,j = c,t incor¬ 
porate QCD corrections of higher order in ag- These 
are obtained by resumming large logarithms with the 
renormalization group evolution m- To make it scale 
and renormalization scheme independent, the renormal¬ 
ization group running factor with 3-flavors 6(/r) is fac¬ 
tored out, 

= (52) 

where (/i) is given in Eq. (A28) of Appendix 

It is combined with the hadronic matrix elements of 
the four fermion operator and used to define 

a renormalization group invariant quantity Bk, 


where 


Bk{h) = 


Bk = BK{H)b{H) 1 

|(ArO|s7^ 75^10) (Olsy^^yscllArO) 


(53) 


(54) 


can be calculated from lattice QCD at a common scale 
such as /r = 2 GeV. Fk is the kaon decay constant. 
Inserting Eq. (48) into Eq. (47), we can identify the 


short distance contribution to m! as follows, 




(55) 


where 


Asd = ImAt 


R.eAc??cc-S'o(a;c) - BeXt'qttSo{xt) 


- (ReAc - ReAt)? 7 ct-S'o(xc,xt) 


(56) 


Here, we use another unitarity identity, ImAj = —ImAc. 
It can be shown from the unitarity condition of Eq. (51) 


and noting that A^ is real in the standard parametriza- 
tion. 

With the Wolfenstein parametrization for the CKM 
matrix elements 1301. 


ReA, = -A ( 1 - y 


l-^-AV(l-p) 


(57) 

(58) 

(59) 

(60) 


They are accurate to 0(A®). Here, we have neglected 
terms of 0{X'^). Then 


ReA* = - (^1 - y ) A2 a5(I - p), 

ImAt = pA^X ^, 
where 


p = p 1- 


r? = r? 1 - 


-^SD = 7A^|I4h|^ 


|f4fep(l - p)'quSo{xt){l + r) 
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+ ( 1 ^ ) {'nctSo{xc,Xt) - riccSo{xc)} 


(61) 


Hence, we take the experimental value of AMk for 


AM in Eq. (46). This is a very good approximation, 


where r = {r]ccSo{Xc) - 2r]atSo{Xc,Xt)}/{r]ttSo{xt)}. 
Here, note that we replace A by Vet, using the relation 
IKhI =41A2 + 0(A®). 


D. Long Distance Contribution 


In the previous section, Section |H C[ we explain the 
short distance contribution of the effective Hamiltonian 
"Hgp to m'. Here, we would like to address the effect of 
the long distance contribution to m'. 

The parts of second order in in Eq. ^ and Eq. © 
correspond to the long distance contributions. The long 
distance contribution TOlq of m' is 


^LD — 


{Ki\H^\C){C\H^\K2) 


= —Im 




c 


nriKO - Ec 

{^\H^\C){C\H^\K^) 
rriKO — Ec 


+ 5mLj3 , 




c 


rriKO - Ec 


AMk = Ml — Ms = Re(AL — As) 


= Re^(AM+^Ar) (l-4e7o)) 
(1 + J cot 0 ) Y^l - 4e^p 


= AM ■ Re 

= AM(l-2Re e(o)(l+ zcot 6 ») + 0 ( 7 ( 0 ))). 


( 66 ) 


Here, note that 9 = ir/A and £(□) = Sk- Hence, the 
difference between AMk and AM is of 0{AMe^Q^). This 
small correction can make a change of O(e^o)) in Ek- 
Here, note that 0 ( 7 ( 0 )) ^ O{oje'). Hence, this is so 
small that we neglect it. 


E. Master Formula: ek 


From Eqs. (46), ( [55| ), (61), and (62), the phenomeno¬ 
logical expression for the indirect CP violation parameter 
in the SM is 


Here, note that vanishes due to CPT invariance. 


drriLY) = 0. 


(63) 


The absorptive part 7 ', which comes entirely from the 
long distance effect, is treated in the previous section. 
Section III HI 

The net contribution ^ld to Ek in Eq. (|^ , which comes 
from m^jo, was estimated to be the same order of mag¬ 
nitude as ^0 using chiral perturbation theory HO]. They 
claim that ^ld = —0.4(3)^o and that ^ld is at most a 
4% correction to Ek- This claim is consistent with the 
estimate of about 2% in Ref. [7]. 

Following the estimate of the long distance contribu¬ 
tion TOlq, it was claimed in Ref. |in] that this contri¬ 
bution should be incorporated. Here, we treat the long 
distance effect of TOlq as a systematic error in the error 
budget of Ek ■ 

The theoretical expression for the mass difference AM 


defined by Eq. (22) is 




AM = 2Re{K''\H^s^\K^>) + AMld 


(64) 


AMld 


2Re 


^ rriKO - Ec 


(65) 


There has been an attempt to calculate AMld in lattice 
QCD [7l|42]. Since the precision of lattice results is not 
as good as that of experiment, we use the experimental 
results for AMk in this paper. 


Ek =e*®-\/2sin0^C'eilif Asd + + ■^ld) 

+ O(w£') + o(?or 2 /ri), 


(67) 


where 


GlFl,mKoM^ 

6V2tt^AMk 


^/2AMK ' 


( 68 ) 

(69) 


Here, ^ld is the long distance effect of « 2% [7]. Precise 
theoretical evaluation of ^ld from lattice QCD is not 
available yet. Hence, we do not include this effect in the 
central value of Ek, but we take it as a systematic error 
in the error budget of Ek- 

The correction terms 0{uje') and Ci(^oP 2 /ri) are of 
order 10 “^, and we also neglect them in this analysis. 

In Eq. (61), the parameter r is very small (« 10“^) 
and also Mj8 « 10“"^. Hence, if we neglect these small 
terms in Eq. (61), we can obtain the same formula as 
in Ref. However, in this paper we keep both the r 
parameter and the A '^/8 term in Eq. (61), even though 
they make no difference to our conclusion. 


In Ref. [S] , the multiplicative factor Kg was introduced 
to incorporate long distance effects ^ld, the small addi¬ 
tive correction and deviation of the angle 9 from the 
value 45°. Since ^0 can be estimated from lattice QCD 
m, we can treat this small contribution to Ek explicitly. 































III. DATA ANALYSIS 
A. Input Parameters 

The CKMfitter and UTfit groups provide the Wolfen- 
stein parameters X,p,fj and A from the global UT fit. 
Here, we use A, p, fj from CKMfitter [55] and UTfit 
[571 155] . and we use Vcb instead of A, Eq. ([M). The 
parameters A, p, and p are summarized in Table D 
The parameters EkjBk, and Vcb are inputs to the 
global UT fit. Hence, the Wolfenstein parameters ex¬ 
tracted from the global UT fit of the CKMfitter and 
UTfit groups contain unwanted dependence on the Ek 
calculated from the master formula, Eq. ( [^ . To self- 
consistently determine Ek-, we take another input set 
from the angle-only fit (AOE) in Ref. [55] . The AOF 
does not use Ek,Bk, and Vcb as inputs to determine the 
UT apex of p and fj |^. The AOF gives the UT apex 
(p, fj) but not A. We can take A independently from the 
CKM matrix element Vus, because this is parametrized 
by 

|K«| = A + 0(A"). (70) 

Here we use the average of results extracted from the 
and 77^2 decays [4]. 


TABLE I. Wolfenstein Parameters 



CKMfitter 

UTfit 

AOF [5^ 

A 

0.22537(61) /[I] 

0.2255(6) /g] 

0.2253(8) /g] 

P 

0.117(21) /[I] 

0.124(24) /g] 

0.139(29) /g3] 

V 

0.353(13) /g] 

0.354(15) /g] 

0.337(16) /gg 


B —>■ D*i9 at zero-recoil {w = 1) [50]. The authors 
of Ref. [50| used the Wilson clover action for the heavy 
quarks, which is tuned by the Fermilab interpretation 
[44] via heavy quark effective theory gSHlZI, with the 
MILC Nf =2-1-1 asqtad gauge ensembles [IS]. The 
heavy quark symmetry and heavy quark effective theory 
play a key role throughout their strategies. Considering 
about a 1% enhancement by the electromagnetic correc¬ 
tion IpemI, they combined their lattice result with the 
HFAG average [49] of experimental values J^(l)|pEM||Kcf)| 
to extract |Ucb|- 


TABLE II. Inclusive and exclusive lUcbl in units of 10“®. 
Here, Kin. represents the kinetic scheme in the heavy quark 
expansion, and IS, the IS scheme. 


Inclusive (Kin.) 

Inclusive (IS) 

Exclusive 

42.21(78) /[H] 

41.96(45)(07) /gD] 39.04(49)(53)(19) /gD] 


There has been significant progress in unquenched 
QCD studies in lattice gauge theory since 2000. This 
progress makes several lattice calculations of Bk avail¬ 
able at A/ = 2-1-1 [l6|[T71[5T|[55|. FLAG provides various 
lattice results for Bk with Nf = 2 + 1 and the lattice av¬ 
erage [T5|. Here, we use the Nf = 2 + 1 FLAG average 
in Ref. [T5| and the SWME result as inputs, which are 
summarized in Table [ml FLAG uses the SWME result 
from Ref. [T5] , which is not much different from the most 
up-to-date value [18] that we use in this analysis. The 
BMW calculation m quotes the smallest error, and it 
dominates the FLAG average. The SWME result [T5] 
quotes a larger error, and its value deviates most from 
the FLAG average. 


The input values that we use for Vcb are summarized 
in Table m The inclusive determination considers the 
following inclusive decays: B —b Xclv and B —)■ 
Moments of lepton energy, hadron masses, and photon 
energy are measured from the relevant decay. Those 
moments are fit to theoretical expressions which are 
obtained by applying the operator product expansion 
(OPE) to the decay amplitude with respect to the strong 
coupling as and inverse heavy quark mass A/mb- There 
are two schemes for the choice of h quark mass mb in the 
heavy quark expansion: the kinetic scheme and the IS 
scheme |1[55]. We use the value obtained using the ki¬ 
netic scheme [55], which has somewhat larger errors and 
also was updated more recently]^ 

The exclusive determination considers the semi- 
leptonic decay of R to D or D*. Here, we use the most 
up-to-date value from the FNAL/MILG lattice calcula¬ 
tion of the form factor F{w) of the semi-leptonic decay 


^ In Ref. |22| . inclusive Ub is determined using the semi-leptonic 
B decays but not radiative B decays. 


TABLE III. Bk 

FLAG SWME 

0.7661(99) /[H] 0.7379(47)(365) /[H] 


The RBG/UKQCD collaboration provides lattice re¬ 
sults for ImA 2 and m- They obtain (defined in 
Eq. (36)) using the relation 



COs{(j)e' — (j^e) ReA 2 /ImA2 
ReAo \ReA2 



(71) 


In this relation, they use the lattice result for ImA 2 and 
take the experimental values for the remaining param¬ 
eters to obtain ^o- In particular, they use the experi¬ 
mental value of Ek as an input parameter to determine 
However, the error is dominated by the experimental 
error of Re(e'/eii'), which is « 14%. In the numera¬ 
tor, cos{4>ei — 4>e) is approximated by 1, because the two 
phases are very close to each other |4], 


= 43.52(5), 


(72) 
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= 42.3(15). (73) 

The final result for in Ref. [IH] is 

Co = -1.63(19)(20) X 10"^ . (74) 

The magnitude of Cld is about 1.6% [7]. We incorpo¬ 
rate the systematic uncertainty in ek due to neglecting 
Cld by treating Cld as a Gaussian distribution about zero 
with a width of 1 . 6 %, 

Cld = (0±1.6)%. (75) 

The factor rju is given at next-to-leading order (NLO) 

in Ref. Other factors rjct and rj^c are given at next- 
to-next-to-leading order (NNLO) in Refs. [53] and [54] . 
respectively. The NNLO values of rjct and rj^c are larger 
than the NLO results in Ref. [S]: 


772^0 = 0.47(4), (76) 

^NNLO ^ 0.496(47), ( 77 ) 

77 NLO = 1.43(23), (78) 

^NNLO ^ p72(27) . (79) 


TABLE V. QCD corrections 


Input 

Value 

Ref. 


^CC 

1.72(27) 

Appendix 

A 

Vtt 

0.5765(65) 

m 



0.496(47) 

m 


TABLE VI. Other Input Parameters 

Input 

Value 

Ref. 

Gp 

1.1663787(6) x 10"® GeV^ 

Ill 

Mw 

80.385(15) GeV 

HI 

mc(mc) 

1.275(25) GeV 

m 


163.3(2.7) GeV 

m 

e 

43.52(5)° 

m 

m^o 

497.614(24) MeV 

m 

AMk 

3.484(6) X 10"^^ MeV 

m 

Fk 

156.2(7) MeV 

m 


B. Error Estimate 


Here, we quote the NNLO result for rjcc from SWME 
in Table EYI which is a major update to our previous 
analysis [5S]. In the case of rjcc, the NNLO correction is 
as large as the NLO correction. Hence, the convergence 
of the perturbative series in rjcc is in question [54] . 

In Ref. [SS], they claim that the error is overestimated 
for the NNLO value of rjcc given in Ref. [M]. Hence, in 
order to check the claim, we follow the renormalization 
group (RG) evolution for rjcc described in Ref. [53] to 
produce the NNLO value of rjcc- The results are summa¬ 
rized in Table IV In this table, note that the results are 
consistent with one another within the systematic errors. 
Here, “SWME” represents our evaluation of rj^c: which 
is essentially identical to that of Ref. [55] . Details of our 
results are explained in Appendix [^ In this paper, we 
use the SWME result for r^cc to obtain ek- 


TABLE IV. Results of tJcc at NNLO. 


collaboration 

Value 

Ref. 


Brod and Gorbahn 

1.86(76) 

m 


Buras and Girrbach 

1.70(21) 

m 


SWME 

1.72(27) 

Appendix 

A 


The input values for 77 ^- that we use in this paper are 
summarized in Table m 

The remaining input parameters are the Fermi con¬ 
stant Gf, W boson mass Mw, quark masses rriq, kaon 
mass 771 jfo, mass difference AMx, and kaon decay con¬ 
stant Fk- These are summarized in Table jVlj 


We use the Monte Garlo method to obtain the expec¬ 
tation value of Ek, 



xp(x)eic(x) 


1 


Ns 

+ O 

i=l 



where x is a sample vector of the input parameters 
that we describe in the previous section. We generate 
Ng = 100 , 000 random sample vectors x that follow the 
multivariate Gaussian probability distribution p(x) with 
covariance matrix Cij = {SxiSxj),6xi = Xi — (xi), 

p(x) =A7exp(- i(x- (x))^C'-i(x- (x))) , (81) 

where Af is the probability density normalization factor. 
The dimension of a sample vector x is d = 18, which is 
the total number of input parameters which appear in 
the master formula for Ek, Eq. We construct the 

covariance matrix by assuming a correlation Cij between 
parameters Xi and Xj 

Cij — CijCJiCTj , ( 1 ^ Cij ^ 1 ) , (^^) 


and using the mean (xi) and error aj of the input pa¬ 
rameter Xi given in Tables [l)[IT HI V VI Eq. (741 , and 


Eq. (75). When the quoted error is asymmetric, we take 
the larger one as a symmetric error. The actual values of 


the correlation matrix Cy are given in Section IV 


In this numerical study, we used the GNU Scientific 
Library (GSL) [58]. Specifically, we used the pseudo ran¬ 
dom number generator ranlxd2 [59] to obtain uniformly 
distributed random numbers. Then we convert them to 
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the multivariate Gaussian distribution using GSL built- 
in functions. 

To find the contribution to the total error from the 
error in each parameter entering the master formula for 
£k, we use the following error propagation method. For 
/ = |£if I, the variance is 

a^ = ([/(x)-/((x))n = ([<5/(x)n, (83) 


where 


N 




1=1 


Sx 


3 ’ 


(84) 


as a linear approximation. Then, the square of the rela¬ 
tive error is obtained by 


N 

E 


9/(x) 


N 

~ E/ ’ 

i,fc=i 


a/(x) 


dxk 


9/(x) 


dx. 


(x> 


^3 

if) 


(SxjSxk) 

9/(x) 


dxL 


^ k 

( 7 ) 


(85) 


where is again the correlation matrix; by definition the 
diagonal components are always ca = 1. This method of 
error propagation is used to cross-check our Monte Carlo 
result. Indeed, errors estimated by these two different 
methods are consistent with each other. And for the 
error budget in Table |IX[ we quote the fractional error 
for the parameter Xi, which is defined as 


a/(x) 


dxi 


(x> 


if) 


{f? 


( 86 ) 


in percent. 


IV. RESULTS 


Let us define as the theoretical evaluation of \£k\ 
obtained using the master formula, Eq. (671. We define 
as the experimental value of \£k\, given in Eq. ([^. 

Let us define Ae^ as the difference between and 

„SM. 

■ 


Ac — e-SM 

L\£k = ~ Sj- 


-K 


(87) 


Here, we assume that the theoretical phase 0 in Eq. (67) 
is equal to the experimental phase 4>e in Eq. § i- 

we present results for e™ obtained using 


In Table 


VII 


the FLAG average for m and Vcb from both inclu¬ 
sive [55] and exclusive channels |5D]. The corresponding 
probability distributions for and are presented 
in Fig. The corresponding results for Aek are pre¬ 
sented in Tab le IVIIIl 


VIIll we find that with inclusive Vcb is 


From Table 

consistent with within Icr. In other words, Ae^ is 


TABLE VII. in 

units of 10 Here. 

we use the FLAG 

average for Bk in Table III The input methods of CKMfitter, 

UTht, and AOF represent different inputs for the Wolfenstein 
parameters, which are explained in detail in Section |IIIA1 

Input Method 

Inclusive Vcb 

Exclusive Vcb 

CKMfitter 

2.31(23) 

1.73(18) 

UTfit 

2.30(24) 

1.73(19) 

AOF 

2.15(23) 

1.61(18) 

TABLE VIII. Ask- 

Here, we use from Table VH We 

obtain a by combining the errors of and in quadra- 

ture. 



Input Method 

Inclusive Vcb 

Exclusive Vcb 

CKMfitter 

-0.34(7 

2.7a 

UTfit 

-0.31(7 

2.7a 

AOF 

0.33(7 

3.4(7 


consistent with zero with inclusive Vcb regardless of the 
input methods. 

However, from Tables 


VH 


and 


VHI 


e™ with exclu¬ 


sive Vcb is only 72% of ■ For this case, with the most 
reliable input method (AOF), Ae^ is 3.4cr. Since the 
largest contribution in our estimate of that we ne¬ 
glect is much less than 2%, the neglected contributions 
cannot explain the gap Aek of 28% with exclusive Vcb- 
Hence, our final results for Ae^^ are 


Aek = 3.4(7 (exclusive Vcb ), (88) 

Aek = 0.33(7 (inclusive Vcb ), (89) 

where we take the AOF result as our quoted value. 

In the case of the FLAG Bk, the BMW result for Bx 
dominates the FLAG result, and the gauge ensem¬ 
bles used for the BMW calculation are independent of 
those used for the determination of exclusive Vcb [2Q] by 
the FNAL/MILG collaboration. Hence, we assume that 
we may neglect the correlation between the FLAG Bk 
and the exclusive Vcb- However, the SWME Bk calcula¬ 
tion in Ref. [T8| shares the same MILG gauge ensembles 
with the exclusive Vet determination in Ref. [5D]. Hence, 
in this case, we cannot neglect the correlation between 
the SWME Bk and the exclusive Vcb- We introduce 
-1-50% correlation and —50% anti-correlation between the 
SWME Bk and the exclusive Vcb-, and take the maximum 
deviation from the uncorrelated case as the systematic er¬ 
ror due to the unknown correlation between them. The 
details of this analysis are explained in Appendix]^ How¬ 
ever, this analysis shows that the size of the ambiguity 
due to the correlation between the SWME Bk and the 
exclusive Vcb is so large that we can use the results of 
the SWME Bk only to cross-check those with the FLAG 
Bk- This analysis of the correlation is another update 
from the previous paper [55] . 














































11 


2a 3a Aa 




1 1.5 2 2.5 

(d)CKMfitter 


2a 3a Aa 




2a 3a Aa 



(c)AOF 



1 1.5 2 2.5 


(f)AOF 


FIG. 2. Gaussian probability distributions for (blue dotted line) and (red solid line). 


using the FLAG Bk- The results of 2(a) 
obtained using the exclusive Vcb- 


2 (b) 


and 


2(c) are obtained using the inclusive Vcb- Those of 2(d) 2(e) 


Here, the results are obtained 
and ! 


2 (f) 


are 


4.0 


3.5 


b 

'5 3.0 

lu 

<1 

2.5 


2.0 


FIG. 3. Recent history of Aex along with the theoretical 
progress. 


NNLO»?cc - 

updated Vcb i 
FLAG Bk \ 

NNLO i]ct \ 


RBC/UKQCD ^0 
LLV Avg. Bk 


2012 2013 2014 2015 

T 


It is interesting to understand the historical evolu¬ 
tion of /S.£k I O' along with the theoretical progress in lat¬ 
tice QCD and perturbative QCD. In Fig. we present 
/S.£kIo as a function of time. In 2012, the RBC/UKQCD 
collaboration reported in Kef. [TS]. In addition to this, 
using the LLV average for Bx [SS], the SWME collab¬ 
oration reported = 2.5cr in Ref. [ST] in 2012. In 


2014, FNAL/MILC reported an updated Vcb in the ex¬ 
clusive channel. Using the FLAG average for Bk [H] 
and the NNLO value of lyct [53|, the SWME collabora¬ 
tion reported the updated = S.Ocr in Ref. |SS] in 
2014j^ In this paper, we investigate issues in the NNLO 
calculation of rjcc O [56] and use the SWME result in 
Table |V| to report the updated Ae^^ = 3.4cr in Eq. (88). 


V. CONCLUSION 

In this paper, we observe that there is a substantial 
3 .4(7 tension in £k between experiment and the SM the¬ 
ory with lattice QCD inputs. For this claim, we choose 
the angle-only fit (AOF), the exclusive Vcb (lattice QCD 
results), and the FLAG Bk (lattice QCD results) to de¬ 
termine the final value. We choose the AOF method to 
determine the final result because the AOF Wolfenstein 
parameters do not have unwanted correlation with £k 
and Bk- However, the tension disappears in the case of 
inclusive V/h (results of the heavy quark expansion based 


^ Here, we evaluate the ek using the correct master formula in 
Eq. ||67^ with the same inputs as in 1611 (2012) and 1551 (2014). 
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TABLE IX. Error budget for obtained using the AOF 
method, the exclusive Vc.b, and the FLAG Bk- Here, the 
values are fractional contributions to the total error obtained 
using the formula in Eq. (|86|. 


source 

error (%) 

memo 

Vcb 

39.3 

FNAL/MILC 

V 

20.4 

AOF 

'fjet 

16.9 

c — t Box 

rice 

7.1 

c — c Box 

P 

5.4 

AOF 

mt 

2.4 


Co 

2.2 

RBC/UKQCD 

Cld 

2.0 

RBC/UKQCD 

Bk 

1.5 

FLAG 

me 

1.0 



on the OPE) regardless of the choices for the Wolfenstein 
parameters_; _ 

we present the error budget of for the 


In Table 


IX 


central value. This is obtained using the error propaga¬ 
tion method explained in Section |III B[ From this error 
budget, we find out that I4b dominates the error in e^. 
Hence, it is essential to reduce the error of I4h as much 
as possible. (See also Refs. [62l[63].) In order to achieve 
this goal, we plan to extract t4b from the exclusive chan¬ 
nel using the Oktay-Kronfeld (OK) action for heavy 
quarks to calculate the form factors for B —)■ de¬ 

cays. Preliminary results in the early stage of the Vcb 
project are reported in Refs. [55H57] . 

Our results for rjcc are consistent with the conclusion 
of Ref. [SI] regarding the convergence of perturbation 
theory. Uncertainty due to truncated higher order terms 
requires further investigation in the future. Lattice QCD 
calculations with dynamical charm quarks, such as that 
envisioned by the RBC/UKQCD collaboration, could 
shed light on this issue. 

We expect that our results for ek would be consistent 
with those from a global UT analysis, such as that in 
Ref. [5D]. The authors of Ref. [SD] performed one analysis 
with inclusive Vcb and another analysis with exclusive Vcb 
instead of inflating their errors and taking the average. 
In this respect, the analysis of Ref. (50] is different from 
those of UTfit and CKMfitter. Such a global analysis 
with up-to-date inputs from lattice QCD has not been 
performed yet. It would be interesting to see the results 
of such an analysis. 
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Appendix A: Next-to-next-to leading order r]c 


We will begin from the master formula for rjcc |54j and 
give an explicit expression for each component which is 
necessary for a numerical evaluation. For fJ. < fic. 


1 


Vcc = 


C-s‘2i^^c) K (Mc,3)]“+(")K;i(/i„3). 


' {me) 


(Al) 

The magic number a+(3) = 2/9 can be obtained from 
Eq. (A32). as(/r,/) is the running strong coupling con¬ 


stant with / active flavors at scale fi. We will use the 
four-loop as running formula HOJ jSS] ■ The Wilson coef¬ 
ficient Cg 2 {t^c) of the AS = 2 four-fermion operator is 
defined by Eq. (jA^ . The running matrix iL_|_^(/ic, 3) is 
given by Eq. ( |A29 ). 

At the charm scale ^c, the effective four flavor theory is 
matched to the effective three flavor theory by requiring 
the following condition m, 

= ^C^s‘2{Qs2) . (A2) 




The matrix elements and the Wilson coefficients are ex¬ 
panded in the three flavor strong coupling as{nc,i)j 


{Qs2) =fis2(Qs2)® , 


fs2 =1 + 


Os (Mo 3) / Os (/^C; 3) 


47r 


' S2 


47r 


=(2) 
' S2 


(A3) 

, (A4) 


(Q^Qj) =- 






d^J{Qs2)^°^ 


(A5) 


d m ^ ( 1 ) ^ / asiiic,^) 

■> 47r \ 47r 


d) 


( 2 ) 


(A6) 


CM =C'f)(Mc) + °%^C'f)(Mc) 

Os(Mei 3) ^ (2) 


47r 


cr{^^c), 


(A7) 


CM) =Clf\t^c) + 


s(Mcj 3) 
47r 


C^fH/re). (A8) 
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Then, the matching results are 

+ (Cf + Cf^cf^)df^ 


(A9) 




(AlO) 


C'sT\dc) =ml{nc) 


(df^ 


oiog^^: 


dc .( 1 ) 


2 U 


+ (cf( d 


( 1 ) 

U 


? log 

3 ^m2 


+ (cf + cf ^)4 


'f(l) ^(1) 
''j 


(0)^(2)>^ 


-j(O) 


(All) 


where rric = mc{mc) in the logarithms multiplied by 


dl°'^\ and 


d(i) =rf(i) _ c7(o)f(i) 


i-j 

J(2) 


j(2) 


^(1)^:(1) _ u(0) ~(2) 


^ij -^ij % ^S2 % ^S2 ■ 


(A12) 

(A13) 


Note that the matching scale is the charm quark mass 
/ic = mc(mc); inEqs. (A9), (AlO), and (All), the Wilson 


coefficients cf\^c) (^ = 0,1, 2; i = ±) are evaluated at 
Me = rric, 




TT 

18 ’ 


Ai) _ 

^S2 — ■ 


17 


(A16) 


(A17) 


The next-to-next-to-leading order (NNLO) calculation 
results are presented in Ref. [54], 


P'’ - 


1665873233 1573 


133 


8164800 


162 72 


, 49^ , , 4313,2 15059, 

36'"^ ~ 1296 

210213,^ 1501^2 7567241^ 

560 ^ ~ ^4“^2 - 204120 

1697893, 11575, 

C3 + ^^C4, (A18) 


7776 


216 


,q2) ^ 87537463 685 83 

+" "+ 1166400 162 72 

695, , 1475,2 57763 

~ ^IfT ~ 1296 


4797^ 791,2 , 366569, 

80 ~ 29160 

57673 4999 

7776 ~ ^IfT'’"^ ’ 


(A19) 


= C)‘^(mc). (A14) 

These are obtained by renormalization group evolution 
from the scale fj,w down to the scale Mc = 'me- (See 
Eq. (A45).) To examine the size of residual scale depen- 


M), 


dence, we vary /ic, keeping the condition Eq. (A14). 

Then the residual scale Mc dependence in Cgl{fj,c) 
enters from logarithms which are shown explicitly in 
Eq. (All) and through df^ and r^ 2 i it also comes from 


the expansion mc{fJ,c)- The expansion of the charm quark 
mass mc{^ic) near = Wc is given by Eq. (A48l with 


/ = 4. The resulting residual scale dependence in rjee can 
be seen from Fig. 

The leading and next-to-leading order (NLO) calcula¬ 
tions can be found from Ref. |69| . with the number of 
colors Nc = 3, lc = log {fil/ml (/ic)), 


- - 

++ “ 2 ’ 

^(0) _^(o) __i 
“+- — «-+ — 2 ’ 

= - 

-- 2 ’ 


(A15) 


-(2) 2129775941 491 „ 11 

dll =-xr::vxrr— + + —D3 


8164800 


162 


72 


865,, 12533,2 171121, 

T -C2^C 4“ -4“ - 

36 ^ 216 1296 

59121^ 517,2 , 9261883, 

560 ~ 204120 

411709 7913 

~ 7776 ~ ■ 

Some constants for the master integrals are m 


(A20) 


CI 2 W 


D 3 =6C3 - ^C4 - 6 
B4=- 4(2 In^ 2 -1 I ln'‘ 2 - ^(4 + 16Li4 f l 


(D 


52 = 


9^3 


CI 2 W 


(i)' 


with 


Cl 2 (:r) =Im (Li 2 (e“)) 

00 


(A21) 

(A22) 

(A23) 


k=l 


++ 2 6 


dE = d'Li = -6Ze - ^ 


W- 


7 ( 1 ) 


23 




57r^ 
~18 ’ 


and the Riemann zeta function is 


Cn = E 


fc=l 


(A24) 






































14 


In numerical evaluation, we use approximated numbers 
which are obtained using Mathematica. 


C 2 =1.644934... = — , 

6 


Ca =1.202056... , 
C 4 =1.082323... = 


90 


Li4 ( i ) =0.5174790... 


Ch 


(i) 


=1.014941... . 


(A25) 


For C 2 and C 4 , we also give the exact expression. 

The renormalization group evolution of the Wilson co¬ 
efficients C± is described by the evolution matrix Uij 

[zn, 


(M? (Mo) 


(A26) 


which is diagonalized by the specific choice of evanescent 
operators, 






(A27) 


where 


^±(m)=i + 


j{2) 


47r 




j: 


± I 


simJ) 


{Jf - 


and 


j(i) ^ 

^ / 3 o ^ 2^0 ’ 


7(2) _ -^2 ,^f( jW\2 / 3 l t(1)\ 


(A28) 

(A29) 

(A30) 

(A31) 


The expansion coefficients of the QCD beta function /3j 


are given in Eq. (A46 1 . The anomalous dimensions 7 ,^ 
for the operators Q± are taken from Ref. [72] . 


i) 


7^°^ = ±6 (1 T ^ ) = 2/?oa± 




‘^5 


7f = ^(349049 ±201485) 


(A32) 

(A33) 

(115577 T 9795)/ 


130 


27 


1 


T 1 ffi F / ffi 672 + 80 1 T 7 / Cs . 


(A34) 


The number of active flavors is fixed while applying 
Eq. (A26). The number of flavors is implied by the strong 


coupling constant in Eq. (A27). 


The initial conditions for the Wilson coefficients C± 
are chosen at the scale 


c±M (m) + 

0‘s{f^W,5)\'^ ^{2) 


477 




(A35) 


The expansion coefficients are given in Ref. 


C7(«) = ±^lT3j (11401.^ 


C^\fiw) = - 


3600 


(135677 T 124095) 






(A36) 


where 


T{xt) =-^ + 32a:t -I- -I- 16a;t^ Inxt 

— (8 -b 16a:t)-\/4xt — ICI 2 ^2 arcsin 


V2v^ 


(A37) 


Xt = and Cl 2 (a;) is given in Eq. (A22). 

In numerical evaluation, we use an approximated num¬ 
ber which is obtained using Mathematica, 


CI 2 ( 2 arcsin 


V2v^ 


= 0.8464504... 


(A38) 


The value Xt is evaluated with the top quark mass 
ixLt{mt) = 163.3 GeV and Mw = 80.385 GeV, approxi¬ 
mating mtifJ-w) = 

Running from to the bottom quark threshold is 
achieved by 

Q(Mfc,5) =K,{^ib,5) 


(A39) 


5) -b 


s(M6,5) 

«s(A^b, 5)^(1) 


477 


Q '(/ib,5) 


«s(Mbi5)^ (2) 


477 


cr(Mb,5). 


The threshold correction at /if, is given by the following. 
Writing 


ai^ibA) =cf (/7b,4) + 


5 (/^b) 4) 
477 


as(Aib,4) (1) 


477 


cr^(Mb,4) 


C'f(/if„4), (A40) 
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then 


=C'«(m6,5), 


(A41) 


where 


3 mr 




= --\n^c2\^^b,5) 


2 

(A42) 


The definition of r) ' and Sr]^ , and their combination, 


and 

the multiplicative factor of , can be found in Ref. m- 
Running from fib to the matching scale fic = mcijric) 
is achieved by 

C,(Mc,4)=/C.(Mc,4) '^\-\fibA)Cr{fibA) 


asificA) 


—^i ^(Mc,4) + 


a«(Mc,4) (1) 


47r 


Q ^(Mc,4) 


5 (/^c; 4) 
47r 


cf^(Mc,4). 


(A43) 


^ 149753 /1078361 6508, , , 

/33 + 3564C3 - ,,,, + ^Cs / 


V 162 27 

/ 50065 6472 \ ,,, 1093 

+ I + —Ca f + ■ 


V 162 


81 


729 


(A46) 


The NNNLO decoupling relation of the strong coupling 
constant at a flavor threshold fi is |40] 


asAJ- 1) 

47r 

asA,f) _ 
Att 


a-sAJ) 

47r 


otsAA) 


Att 


\\n^ 


22 38 n" 4 2 

v;-In —T + X In —v 


TOt 


ml 


asAA) \ / 564731 
47r y ^ 1944 


2633,, 82043^ 

-\I - - 


486 


432 


+ ;^ln^(-6793 + 281(/-l)) 


131 2 A 

— In — 
9 m 


27 m 


(A47) 


where mh = mh(mh) is the scale invariant MS mass of 
the heavy flavor which is removed from an effective the¬ 
ory below the threshold fi, and ^3 is a Riemann zeta 
function, Eq. (A24|. 


The running quark mass mq{fi), an MS mass at scale 
fi, for a fixed number of active flavors / is |40] 


mg{fi) R{as{fi)/ATr) 


mgifio) R{as{fio)/ATr) 


(A48) 


In the matching calculation, we need the expansion in 
the three flavor strong coupling, Eq. (A71. 


E quatin g Ci(fic) in Eq. (A7) to Ci{fic, 4) after applying 
Eq. ( A47|) to the flavor threshold with f = A, then 


= 6 -® (Me, 4), 


Cf (Me) (Mc, 4) + |ln5cf)(Mc,4). (A44) 

O /7i_ 


(1). 


Hence, at the matching scale of the charm quark mass 
Me = mcirric), we obtain 


C'f(m,)=Cr(me,4), (^ = 0,1,2). 


^(0 


(A45) 


The QCD beta function expansion coefficients j3i to 
four-loop order are gniiig: 


/^o=ll-3/, 

A =102-^/, 


with 


R{x) =x‘^° < 1 + (ci — biCo)x 


1 


+ 2 “ bicof +C 2 - bici + 61 Co - &2C0] x'^ 


+ 


-(ci - biCof 

+ ^(Cl - 5 iCo)(C2 - 61 C 1 -I- 61 Co - 62 C 0 ) 
+ ^(C 3 - biC2 + bjci - 62C1 - &?Co 


26162C0 - 63C0) 


(A49) 


where mq{fio) is the scale invariant mass mg = mq{mq), 

bi = Pi/Po, and Ci = 'ym/Po- The QCD beta function co¬ 
efficients Pi are given in Eq. (|A46 1 . The mass anomalous 


dimensions are known up to four-loop order, 

=4 

(1) 202 20 ^ 

Tm ^ 9 ’ 
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FIG. 4. Scale /Xc dependence of rjcc- Dotted, solid, and dashed 
lines represent results for = 40, 80, and 160 GeV, respec¬ 
tively. The upper three (red) lines are the NNLO results, and 
the lower three (purple) lines are the NLO results. 


7 ^') =1249 - 


/2216 160 \ 1402 


^(3) = 

/ m 


162 


/5242 


— -f - 

27 

91723 

34192 

27 

9 ‘ 

800 . 



C3 - 88OOC5 

Cs + 88OC4 ■ 


18400 


Cs / 


r+ - 


332 64 


Ca 


f- 

(A50) 


The value of is obtained by summing the first two 

terms in the series, and the value of is obtained 

by summing all three terms in the series. 

At the scales = 80 GeV and /Tc = 1.300 GeV, 


Vcc/[as{Hc, 3)]“+ =1.113769 -k 0.568911 -k 0.433783, 

^NNLO = 1 . 720690 , 


r;NLO = 1 . 368023 . 

(A53) 

We claim the NNLO rfee is 


jyNNLO ^ 1.72(27). 

(A54) 


The central value corresponds to the result with the 
scales = 1-3 GeV and = 80 GeV. For the er¬ 
ror, we add the /ic and /rvv dependences in quadrature, 

=r;NNL°(/rc = 1-3 GeV,Atw = 80 GeV) 

- = 1-8 GeV, fiw = 80 GeV) 
=0.266, 

^tJ.w (fj'c = 1-3 GeV,fiw = 160 GeV) 

- = 1.3 GeV, nw = 40 GeV))/2 

=0.047. (A55) 


Though we include errors from the in puts as{Mz, 5) and 
rricirnc), the final errors in Eq. (A54) are not affected; we 
use the errors quoted in Ref. 


Sc, = 0.06, S„,, = 0.01. (A56) 


In numerical evaluation, we use approximated numbers 
for the Riemann zeta functions Cn which are obtained 
using Mathematica, 


In Ref. |54j, the authors also added the absolute shift 
from the NLO value of r]cc- It is the main reason for their 
larger error. 


Cs =1-036927... , (A51) 


,NNLO ^ 1,87(76). (A57) 


and ^3 and (4 are given in Eq. (A25). 

We used Eq. (A48) to expand the charm quark mass 
about rric = mc{mc) with / = 4. 

Here, we will give numerical results for an initial 
scale /iw = 80 GeV and a varying charm scale ^c, 
1.0 < He < 2.0 GeV. To examine the dependence 
on the scale Hw^ we repeat the same analysis with 
Hw = 40,160 GeV. The dependence on Hb and mt(mt) is 
ignored [54] . Fig. [^summarizes these results. The follow¬ 
ing are kept fixed for all analyses: the gauge boson masses 
Mz = 91.1876 GeV ,Mw = 80.385 GeV; quark masses 
= 163.3 GeV ,mb{mb) = 4.163 GeV ,mc{mc) = 
1.279 GeV; bottom quark threshold fib = 5-0 GeV; and 
the strong coupling constant that provides an initial value 
for the running formula, as{Mz, 5) = 0.1184 GeV. 

At the scales fiw = 80 GeV and fic = 1.279 GeV, 


The main concern for adding this shift is poor conver¬ 
gence of the as expansion for rfee- We, however, take 
the view that the error from fic dependence suffices to 
estimate the size of higher order corrections. 

Here, we would like to comment on our choice of 
1.3 GeV < He ^ 1.8 GeV. If we examine each order in 
the perturbative corrections of r]cc, 

Vec = l(LO) + 0.37 (nlo) + 0-36(nnlo) 

= 1.72 (SWME), (A58) 

Vec = l(LO) + 0.38 (nlo) + 0.49(nnlo) 

= 1.87 (Brod & Gorbahn). (A59) 

The error that we quote is 0.27, which is roughly equiva¬ 
lent to the size of the NNLO correction^ If we take the 


Vcc/[as{Hc, 3)]“+ =1.129757 -I- 0.571608 -f 0.430890, 

^NNLO = 1 . 738396 , 

^1.387098. (A52) 


^ If we use the size of the NLO correction given in Ref. [^, the 
NNLO correction is 0.29, which is roughly equivalent to our error 
estimate 5^ 0.27. 
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interval 1 GeV < < 2 GeV as suggested in Ref. [M] , 

the error becomes 0.42, which is larger than the NNLO 
correction. If we follow the procedure of Ref. [^, we add 
the NNLO size of 0.36 to this in quadrature which leads 
to a total error of 0.55. This value is significantly larger 
than the size of the NNLO correction. This estimate be¬ 
comes even larger than that of the NLO and NNLO cor¬ 
rections combined in quadrature. If we assume that per¬ 
turbation theory is working, then it is already arguably 
conservative to choose the size of the NNLO correction 
as the systematic error due to truncated NNNLO terms. 
Hence, we believe that the error quoted in Ref. [M] is 
somewhat overestimated. 

In addition, in Ref. [SS] , Buras and Girrbach suggested 
that, if one chooses the interval 1.3 GeV < < 1-8 GeV, 

one can obtain their result for 

« 1.70(21), (AGO) 

which is obtained indirectly through an estimate of the 
long distance contribution to AMk based on a large N 
QGD inspired model calculation. We have directly veri¬ 
fied their claim. 

Our result for agrees with the value quoted in 

Ref. [M] 

jyNLO ^ i.38(52)(07)(02) . (A61) 

Appendix B: with the SWME Bk 


ing results for Aex- In Fig. we show the probability 
distribution for the corresponding e|^. 


TABLE X. in units of 10"^ We use the SWME Bk 
with no correlation between Bk and Vcb- 


Input Method 

Inclusive Vcb 

Exclusive Vcb 

CKMfitter 

2.22(25) 

1.66(19) 

UTfit 

2.21(25) 

1 .66(20) 

AOF 

2.07(25) 

1.55(19) 


TABLE XI. Aek with 
We take from Table 

no correlation between Bk and Vcb- 
is given in Eq. Q . 

Input Method 

Inclusive Vcb 

Exclusive Vcb 

CKMfitter 

0.04a 

2.9a 

UTfit 

0-06a 

2.9a 

AOF 

0.65cr 

3.5a 


TABLE XII. in units of 10"F We use the SWME Bk 
and the exclusive Vcb with -1-50% correlation and —50% anti¬ 
correlation between them. 


Input Method 

c = -50% 

c = -j-50% 

CKMfitter 

1.66(17) 

1.67(22) 

UTfit 

1.66(17) 

1 .66(22) 

AOF 

1.55(17) 

1.55(21) 


Lattice results for the exclusive Vcb 110] and the SWME 
Bk [H] are obtained using overlapping subsets of the 
MILC asqtad gauge ensembles [48]. This implies that 
there exists a complicated, non-trivial correlation be¬ 
tween them. It is possible to calculate, in principle, this 
correlation exactly from the data set. Unfortunately, this 
correlation is not available yet. Hence, the current situa¬ 
tion is that we need to estimate the systematic error due 
to the unknown correlation between Bk and Vcb- 

Here is our strategy. First, we take the uncorrelated 
case as the central value. Second, we introduce -1-50% 
correlation between Bk and Vcb and obtain results for 
Third, we introduce —50% anti-correlation be¬ 


-SM 

'-K 


tween Bk and Vcb and repeat the analysis to obtain e^. 
Fourth, we take the maximum deviation from the cen¬ 
tral value as the systematic error due to the unknown 
correlation between Bk and Vcb- 

In Table we present results for for the uncor¬ 
related case. In Table jXll w e present the corresponding 
results for Aek- In Fig. we show the corresponding 
probability distr ibution of . 

In Table XH we present results for e|U with -1-50% 
correlation and —50% anti-correlation between Bk and 
(exclusive) Vcb- In Table XHI we present the correspond- 


TABLE XIII. Aek with -1-50% correlation and —50% anti¬ 
correlation between Bk and exclusive Vcb- We take e|U from 
Table XII is given in Eq. (|^. 


Input Method 

c = -50% 

c = +507 o 

CKMfitter 

3.4cr 

2 .6a 

UTfit 

3.3cr 

2.5a 

AOF 

4.1a 

3.1a 


Hence, we obtain the final results for Ae^ for the 
SWME Bk and exclusive Vcb- 


Ask = (3.5±0.6)ct, (HI) 


where the error represents the uncertainty due to the 
correlation between Bk and exclusive Vcb- 

First, the results in Eq. (B11 are consistent with those 
in Eq. (88) within the systematic errors. Second, the 
correlation between Bk and exclusive Vcb dominates the 
error in Aek with the SWME Bk- In addition, this error 
is much larger than that in our final results in Eq. ( 881 . 


Hence, we use the results with the SWME Bk only to 
cross-check those with the FLAG 13k- 
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2c7 3(t 4a 



2a 3a 4a 



1.5 


2 2.5 

(b)UTfit 




1.5 2 

(e)UTfit 


2.5 



1.5 


2 2.5 

(c)AOF 



1 1.5 2 2.5 

(f)AOF 


(d)CKMfitter 

FIG. 5. Gaussian probability distributions for (blue dotted line) and (red solid line) with the SWME Bk- Here, we 


assume no correlation between Bk and Vat- Results of 5(a) 5(b) and 5(c) are obtained using the inclusive Vcb- Results of 5(d) 
5(e) and|5(f)|are obtained using the exclusive Vcb- 


[1] J. Ghristenson, J. Cronin, V. Fitch, and R. Turlay, 
Phys.Rev.Lett. 13, 138 (1964). 

[2] A. Alavi-Harati et al. (KTeV Collaboration), 
Phys.Rev.Lett. 83, 22 (1999), hep-ex/9905060. 

[3] V. Fanti et al. (NA48 Collaboration), Phys.Lett. B465, 
335 (1999), hep-ex/9909022. 

[41 K. Olive et al. (Particle Data Group), Chin.Phys. C38, 
090001 (2014). 

[5] M. Kobayashi and T. Maskawa, Prog.Theor.Phys. 49, 
652 (1973). 

[6] N. Cabibbo, Phys.Rev.Lett. 10, 531 (1963). 

[7] N. Christ, T. Izubuchi, C. Sachrajda, A. Soni, and J. Yu 
(RBC and UKQCD Collaborations), Phys.Rev. D88, 
014508 (2013), 1212.5931. 

[8] N. Christ, T. Izubuchi, C. T. Sachrajda, A. Soni, 
and J. Yu (RBC, UKQCD), PoS LATTICE2013, 397 
(2014), 1402.2577. 

[9] A. J. Buras and D. Guadagnoli, Phys.Rev. D78, 033005 
(2008), 0805.3887. 

[10] A. J. Buras, D. Guadagnoli, and G. Isidori, Phys.Lett. 
B688, 309 (2010), 1002.3612. 

[11] J. Donoghue, E. Golowich, and B. R. Hol¬ 
stein, Dynamics of the standard model 
(Camb.Monogr.Part.Phys.Nucl.Phys.CosmoL, 1992). 

[12] S. Aoki, Y. Aoki, C. Bernard, T. Blum, G. Colangelo, 
et al. (2013), 1310.8555. 


[13] T. Bae, Y.-C. Jang, H. Jeong, J. Kim, J. Kim, et al., PoS 
LATTICE2013, 476 (2014), 1310.7319. 

[14] R. Arthur et al. (RBC Collaboration, UKQCD Collabo¬ 
ration), Phys.Rev. D8T, 094514 (2013), 1208.4412. 

[15] J. Laiho and R. S. Van de Water, PoS LATTICE2011, 
293 (2011), 1112.4861. 

[16] T. Bae et al., Phys.Rev.Lett. 109, 041601 (2012), 
1111.5698. 

]17] S. Durr, Z. Fodor, C. Hoelbling, et al., Phys.Lett. B705, 
477 (2011), 1106.3230. 

[18] T. Bae et al. (SWME Collaboration), Phys.Rev. D89, 
074504 (2014), 1402.0048. 

[19] T. Blum, P. Boyle, N. Christ, N. Garron, E. Goode, et al., 
Phys.Rev.Lett. 108, 141601 (2012), 1111.1699. 

[20] J. A. Bailey, A. Bazavov, C. Bernard, et al., Phys.Rev. 
D89, 114504 (2014), 1403.0635. 

[21] P. Gambino and C. Schwanda, Phys.Rev. D89, 014022 
(2014), 1307.4551. 

[22] A. Alberti, P. Gambino, K. J. Healey, and S. Nandi, 
Phys.Rev.Lett. 114, 061802 (2015), 1411.6560. 

[23] N. Uraltsev (2000), hep-ph/0010328. 

[24] G. Buchalla, A. J. Buras, and M. E. Lautenbacher, 
Rev.Mod.Phys. 68, 1125 (1996), hep-ph/9512380. 

[25] J. Charles et al. (CKMfitter Group), Eur.Phys.J. C41, 1 
(2005), hep-ph/0406184. 

[26] A. Hocker, H. Lacker, S. Laplace, and F. Le Diberder, 


































































































































19 


2(7 3(7 4(7 


2(7 3(7 4(7 


jbl. 


1 1.5 2 

(a)CKMfitter 


2.5 



1 1.5 2 2.5 

(d)CKMfitter 





1.5 2 

(c)AOF 


2.5 


2(7 3(7 4(7 



FIG . 6. Gaussian probability distributions for and with the SWME Bk and ex clusi ve Vc b- Results of 6(a) 6(b) 
and 6(c) are obtained with —50% anti-correlation between Bk and Vet- Those of 6(d) 6(e) and 6(f) are obtained with -1-50% 
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